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Abstract 

It is shown that a sufficiently strong four-spin interaction in the spin- 1/2 
spin ladder can cause dimerization. Such interaction can be generated either by 
\ phonons or (in the doped state) by the conventional Coulomb repulsion between 

q | the holes. The dimerized phases are thermodynamically undistinguishable from 

the Haldane phase, but have dramatically different correlation functions: the 
dynamical magnetic susceptibility, instead of displaying a sharp single magnon 
peak near q = ir, shows only a two-particle threshold separated from the ground 
state by a gap. 

PACS numbers: 75.10 Jm, 75.40 Gb. 

The qualitative difference between the universal properties of one-dimensional 
Heisenberg antiferromagnets with half-integer and integer spin, originally predicted 
by Haldane is now well understood. Typical examples are the S = 1/2 and S 
= 1 chains. The spin S = 1/2 chain has a quasi-ordered singlet ground state with 
algebraically decaying spin correlations. Its spectrum is gapless, and the elementary 
excitations (spinons), carrying spin 1/2, in all physical states with integer total spin 
appear only in pairs @]. Deconfinement of the spinons implies that conventional S 
= 1 magnons fail to be stable quasiparticles: the spectral density of the staggered 
magnetization n(x) ~ (— l) n S n shows a purely incoherent background. 
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An alternative picture of a disordered spin liquid with a Haldane gap is effectively 
realized when two S = 1/2 Heisenberg chains are put together to form a spin ladder 
0. In the standard model of spin ladders the interchain interaction is the Heisenberg 
exchange Jj_; at J_l ^ the spinons confine to form triplet (magnon) and singlet 
excitations with gaps m t and m s (mt,m s ~ J±). The triplet excitations contribute a 
coherent 5-peak to the dynamical spin susceptibility x"(Qi u ) near q = it and uj = m t 
HI. In this respect the spin liquid behaves like a conventional magnet with the only 
difference that here the magnons are "optical", that is have a spectral gap. 

In this letter we discuss an example of a disordered spin liquid with a gapful spec- 
trum but without coherent magnon excitations. Such a state can be realized in an 
extended model of the spin S = 1/2 Heisenberg ladder which, apart from the direct 
transverse exchange J±, also includes a four-spin interaction 

w = J \\ E E s »- S i("+ 1 ) 

j=l,2 n 

+ EfJiSxCn) • S 2 (n) + V(S l (n) • S x (n + l))(S 2 (n) • S 2 (n + 1))]. (1) 

It will be assumed that |Jj_|, |V| <C Jn. The new term in (|l|) represents an interchain 
coupling of the spin-dimerization fields e 3 -(n) = (— l) n Sj(n) ■ Sj(n + 1) (j = 1, 2) which 
can be either (i) effectively mediated by spin-phonon interaction (see below) or (ii) in 
the doped phase generated by the conventional Coulomb repulsion between the holes 
moving in the spin correlated background (0). 

In a single S = 1/2 Heisenberg chain described in the continuum limit in 
terms of a massless scalar field $ s (x) ||, the dimerization operator e(x) ~ 
cos ^/2tt& s (x) is a strongly fluctuating field with the same scaling dimension 1/2 as 
that of the three components of the staggered magnetization (— l) n S n ~ n(x) ~ 
(cos \/27iO s (x), sin \^2ttQ s (x), — sin y2ir$ s (x)), where Q s is the field dual to $ s . There- 
fore the interchain coupling term Ve\€2, being as relevant as the transverse exchange 
J^ni • n 2 , may significantly affect the low-energy properties of spin ladders. Indeed, 
while for small \V\ the generalized model (|Tj) occurs in the same Haldane spin liquid 
phase as that of the conventional spin ladder (V = 0), at large enough \ V\ the model ([j]) 
displays non-Haldane spin liquid phases characterized by a spontaneous dimerization. 
The thermodynamic properties of these phases, determined by the absolute values of 
the gaps, are indistiguishable from those for the Haldane phase. However, the cor- 
relation functions of the two systems differ drastically: in a non-Haldane spin liquid 
we are going to discuss the spectral function of n is entirely exhausted by an incoher- 
ent background. This picture should be opposed to the dimerized (spin-Peierls) S = 
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1/2 Heisenberg chain with alternating exchange J n , n +i = J\\ + <K _ l) n ] where, due to 
the presence of the kink-antikink bound states in the spin excitation spectrum, triplet 
magnons constitute well defined, coherent quasiparticles. 

It has been shown earlier || that, in the continuum limit, the Hamiltonian ([[]) 
decouples into four massive real fermionic fields, or equivalently, four noncritical 2D 
Ising models Q, exhibiting the underlying SU(2) x Z 2 symmetry: H = J2 a =i,2,s H mt [£ ] + 
#m s [£o], where 

71) 

H m [£] = —f(.tR d *£* - h d *&) ~ im ^- (2) 
The triplet and singlet masses are given as follows: 

m t = ira (J ± \a s \ 2 -V\a c \ 2 ), m s = -na (3J±\a s \ 2 + V\a c \ 2 ). (3) 

where a c and a s are nonuniversal constants. The thermodynamic properties are deter- 
mined only by \m t \, \m s \, but the symmetry of the ground state and the behavior of the 
dynamical sisceptibility x(li u ) crucially depend on the relative signs of the masses mt 
and m s as well. This follows from the fact that the representation for the staggered spin 
density and dimerization operators in terms of the Ising order and disorder variables 0] 
depends on the sign of mtm s . It will be assumed below that V < 0, while the sign 
of J± can be arbitrary. We shall show that for different signs of mt and m s , as long 
as the triplet branch of the spectrum remains the lowest (\m t \ < \m s \), the two-chain 
ladder is in the Haldane-liquid phase with short-ranged correlations of the staggered 
magnetization and dimerization field, but with coherent S = 1 and S = single-magnon 
excitations. If m t and m s have the same sign, the ground state represents a dimerized, 
spin-disordered phase in which coherent magnon modes are absent. 

Transitions from the Haldane to dimerized phases take place when either the triplet 
excitations become gapless, with the singlet mode still having a finite gap, or vice 
versa. (Notice that in the conventional spin ladder with V — 0, both masses vanish 
simultaneously at J± — ► 0, and the ladder trivially decouples into a pair of independent 
S = 1/2 Heisenberg chains.) The transition at m t = belongs to the universality 



class of the critical, exactly integrable, S = 1 spin chain [ 10 1 ; the corresponding non- 



Haldane phase with \m t \ < \m s \ represents the dimerized state of the S = 1 chain with 



spontaneously broken translational symmetry and doubly degenerate ground state |Tl 

The critical point m s = is of the Ising type; it is associated with a transition to another 

dimerized phase (\m t \ > \m s \), not related to the S = 1 chain. 

1 This decoupling is valid up to a weak residual (marginal) interaction whose role in the gapped 
phases is exhausted by the mass and velocity renormalization [0. 
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We start our discussion by considering the standard J\\-J± two-chain ladder, with 
spin-phonon coupling included via a substitution J\\ — > Jy(n, n + 1) = J» + X[uj(n) — 
Uj(n + 1)]. Only the staggered part of the lattice displacement field along the chains, 
Uj(n), couples to the spin-dimerization operator £j(n). We assume that u ^> Jj_, 
where uq is the phonon frequency at 2k p, as is the case for the most ladder systems 
known. 

In this limit, the 2fcF-phonons can be treated in terms of a quantum Gaussian 
random field whose main effect is mediation of an instantaneous effective coupling 
between the spin-dimerization fields of each chain. We replace the staggered parts of 
Uj(n) by real scalar fields Aj(x), Uj(n) = (— l) n (a/2\)Aj(x), and ignore the kinetic 
energy of vibrations. Then integrating over the displacement fields Aj we get the 
effective interaction 

AS eff [e] = - 9l{ \\ 1 ^ ) J dxdr(el + 4) - / dxdr e x e 2 . (4) 



where go = \/JK\\a is the spin-phonon coupling constant, 7 = K±/K\\ is the ratio of 
the longitudinal and transverse spring constants. The interchain coupling term ~ t\t2 
in (Q) is a relevant perturbation with dimension 1, as opposed to terms ~ e 2 which are 
only marginal and therefore can be neglectedF]. This explains the origin of the extra 
V-term in the extended spin ladder model ([]]) and fixes the value of the constant 

V~-S 2 7 /(2 + 7)<0. (5) 

Let us consider various phases of the model (|l|). We assume that V < 0, while the 
sign of J± can be arbitrary. If Jj_ < 0, the mass m s > 0, while the sign of m t depends 
on the strength of |V| such that for small \V\ the signs of the two masses are opposite. 
All results obtained in Ref.4 are applicable to this case, and we present them briefly 
for completeness. The total (n + = m + 112) and relative (n~ = ni — 112) staggered 
magnetization are given by[|J: 

n + ~ (01//2C3CQ, o-iO- 2 li3(T ), (6) 

n~ ~ (/ii(T 2 /i 3 /i , criyU 2 /i3/^o, ^i^3^o), (7) 

where a a and \i a (a = 1,2,3) are order and disorder parameters of three, SU(2) de- 
generate noncritical Ising models corresponding to the massive triplet of the Majorana 



2 Since the dimcnsionlcss spin-phonon coupling constant <7q jv s is assumed to be small, these terms 
only slightly renormalize the amplitude of the already existing marginally irrelevant perturbation, 
gjft ■ Jl (g ~ 1), to the Gaussian fixed point of the critical S = 1/2 spin chain M. 
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fields £ a , while a , fi Q refer to the fourth, "singlet" Ising model (£ )- The Ising repre- 
sentation for total and relative dimerization fields, e± = e x ± e 2 , can be similarly found 
to be 

6 + ~ yUl/^/^O, C_ ~ a^^flQ. (8) 

The triplet mass determines the deviation of the Ising models from criticality: mt ~ 
(T — T c )/T c < 0; this corresponds to the ordered Ising phase with (/ij) = 0, (<jj) ^ 
(j = 0, 1, 2, 3), in which the two-point correlation functions are given by |8j 



(o-(r)o-(O)) = G ff (mr) ~ A 

A 



Wmr) 2 



2\m\r _|_ /") f „~4|m|r ^ 



(9) 



<Mr)/z(0)> = G At (mr)~-K (|m|r) + O(e- 3 H'-), (1 ) 

7T 

where A is a nonuniversal parameter, i^oO 3 ?) is the MacDonald function, and r = 
(x,v s t). 

In the region of parameters where \m t \ <C |m s |, the lowest (triplet) part of the 
spin ladder spectrum describes universal properties of the S = 1 spin chain with the 
conventional and biquadratic exchange, H = JJ2n[Sn • S n+ i — (3(S n ■ S n+ i) 2 ], near the 
critical (exactly integrable) point (3 = lflKj. This correspondence in valid for any sign 



of m t . The present case m t < describes the Haldane massive phase (/3 < 1), and 
the magnitude of the Haldane gap is given by \m t \ ~ 1-/3 All results obtained 
in Ref.4 are applicable here: the leading asymptotics for the spin correlation functions 
obtained from (H) - (|T0|) 

(n+(r)n + (0)) ~K (|m t |r)[l + O(e- 2|mt|r )], (n-(r)rr(O)) ~ l_ e -W™ t \+m s )r (n) 

reveal the role of n + which determines the staggered magnetization of the effective 
S = 1 chain. Since -Ko(|^l r ) is the real-space propagator of a free massive boson, 
(n + (r)n + (0)) contributes a 5-peak to the imaginary part of the dynamical spin suscep- 
tibility corresponding to a massive triplet magnon. At q ~ 0, x"( w ? Q) is determined 
by the correlations of slow components of the total (J+ = Ji + J 2 = — (z/2)(£ x £)) 
and relative (J_ = Ji — J 2 = iC.C.o) magnetization, giving rise to thresholds at 2\m t \ 
and \m t \ + m s , respectively. Similarly, one finds: 

(e_(r)e_(0)> ~ K (\m s \r) [l + 0(e- 2 "^)] , Mr)e+(0)> ~ -^e" 3 ^, (12) 

implying that the singlet magnon is also a coherently propagating particle, - an elemen- 
tary excitation of the relative dimerization. Excitations of the total dimerization have 
a 3\m t \ threshold. On increasing \V\, \m t \ decreases while m s increases, the inequality 
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\m t \ < m s thus becoming stronger which makes the low-energy effective picture of the 
gapful Haldane phase of the S = 1 spin chain only better. At m t = 0, m s = 4\J±\ the 
triplet of the Majorana fields becomes massless, while the singlet Majorana fermion re- 
mains massive. At this point the correlation functions of the relative staggered magne- 
tization and dimerization field follow power laws: (n - (r)n~(0) ~ (e_(r)e_(0)) ~ r~ 3 / 4 . 
This critical point belongs to the universality class of the level k = 2 SU(2)-symmetric 
Wess-Zumino-Novikov-Witten model with the central charge c = 3/2 ||10|| . 

Further increase of \V\ makes m t positive. The change of sign of m t m s amounts to 
the duality transformation in the singlet (£o) Ising system, implying that in formulas (|6]) 
- (H) the order (ctq) and disorder (// ) parameters must be interchanged. Moreover, since 
we are now in the disordered Ising phase {m t ~ T — T c > 0), < fio >^ 0, < <r >= 0. 
the right-hand sides of formulas (^) and ([H]) should also be interchanged. As a result, 
the spin and dimerization correlation functions are now given by different expressions: 

(n+(r) ■ n+(0)> ~ K 2 (m t r), ( n -(r) • ir(0)> ~ K (m t r)K (m s r) (13) 



(e + (r) e+ (0))~C 



l + 0( 



e -2m t , a 



(e_(r)e_(0)> ~ K^m t r)K (m s r), (14) 



where C is a constant. 

^From ( Il[ ) we conclude that the new phase is characterized by long-range dimer- 
ization ordering along each chain, with zero relative phase: (ei) = (e 2 ) = (l/2)(e + ). 
In the decoupling limit, J± = V = 0, each Heisenberg chain possesses a ^-symmetry: 
this is the symmetry with respect to the interchange of even and odd sublattices gen- 
erated by one lattice spacing translation. The interchain coupling lowers this Z 2 x Z 2 
down to Z 2 , the symmetry under simultaneous translations by a on the both chains. 
Thus, the onset of dimerization is associated with spontaneous breakdown of the resid- 
ual translational ^-symmetry, taking place when the interaction \ V\ exceeds a critical 
value \V C \ ~ Jj_. This phase coincides with the dimerized phase of the generalized S = 



1 spin chain with a biquadratic exchange! 11 



In the dimerized phase the spin correlations undergo dramatic changes. Using long- 
distance asymptotics of the MacDonald functions, from (O) we obtain: 



3m x+(^,vr - q) 



(oj 2 — q 2 — 4m^ 



m t \lu) 2 — q 2 — Am 2 



^mx-(^vr-g) ~ (15) 

y/m t m s ^u z - q z - (m t + m s ) z 

where ±-signs refer to the case where the wave vector in the direction perpendicular to 
the chains is equal to and ir respectively. We observe the disappearance of coherent 
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magnon poles in the dimerized spin fluid; instead we find two-magnon thresholds at 
u = 2m t and to = m t + m s , similar to the structure of x?'(u), q) at small wave vectors 
in the Haldane fluid phase. The fact that two massive magnons, each with momentum 
q ~ 7i, combine to form a two-particle threshold, still at q ~ ix rather than 2n = is 
related to the fact that, in the dimerized phase with 2ao-periodicity, the new Umklapp 
is just 71". 

To get a better understanding of the fact that in the dimerized phase the spectral 
weight of the spin excitations is entirely incoherent, it is instructive to consider the 
limiting case J± = 0. Since €j (j = 1, 2) are scalars in spin space, the model ([[J) displays 
the SU(2) x SU{2) rs SO (A) symmetry with respect to independent spin rotations on 
each chain. Thus, the spectrum is described by 0(4) quadruplet of massive Majorana 
fermions (with the mass m = — |^|(7rao)|a c | 2 ), or equivalently, two noninteracting 
massive Dirac fermions, - quantum solitons of two decoupled sine-Gordon models with 
the coupling constant (3 2 = 4tt: 



The massive Dirac fermions are in fact quantum domain walls (kinks) connecting two 
^-degenerate dimerized vacua of the two-chain system with (e+) = ±|eo|- This can be 
shown by the following simple consideration. A dimerization kink assumes local changes 
6i — > — £i, £2 — > —€2- This is equivalent to simultaneous translations by ao on each 



chain, under which $j — > $j ± y ir/2. This, in turn, reduces to one of two possibilities: 
$+ — > $+ + 0F, $_ —> or $ + — >• $ + , $_ — >• $_ + 0F. But ^Jtx is just the period 
of the cosine potentials in (jl6|). Therefore, a single kink of the relative dimerization is 
nothing but a quantum sine-Gordon soliton (Dirac fermion) carrying a unit topological 
charge, either in the (+) or (— ) channel. Now, in all physical excitations defined in the 
sector with zero total topological charge, the solitons can appear only in pairs. Since 
two or more massive solitons (Dirac fermions) cannot propagate coherently, and since 
there are no soliton-antisoliton bound states at f3 2 = 4n (free fermions!), one has to 
conclude that there will be no particle-like 5-function peaks in the spectral function of 
any local physical quantity of the system. 

Now let us consider the case J± > 0, when m t > while m s may change its sign. 
If m s < 0, the definitions (0)-(f|) still hold, but since m t ~ T — T c > 0, we are in the 
disordered Ising phase. This case can be mapped onto the J± < one by interchanging 
all order and disorder parameters which is equivalent to the replacements n + <-> n~, 
e + <-> e~. As long as \m s \ > m t , the triplet branch of the spectrum describes the 
Haldane phase with the relative magnetization n~, forming effectively the staggered 




7ra 



m 




s 



(16) 
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component of the spin density in the S = 1 spin chain: (n~(r) ■ n~(0)) ~ K (\mt\r). 

On increasing \V\ \m s \ decreases, and the inequality \m s \ > m t will eventually be 
replaced by the opposite one, \m s \ < m t . One might argue that as long as \x\ » £ s 
(the maximal correlation length £ s ~ v s /\m s \), the asymptotic r— dependence of the 
n~-correlator remains intact and the coherent magnon peak, though with a reduced 
amplitude, should still exist. However, this is not so because the large- distance (|r| ^> 
£ s ) asymptotics of the n~-correlator determines its Fourier transform (n~ ■ \i~) q ^ at 
| -7r — q\, \u\ <C m s , so that energies \u\ ~ m t are not accessible. At these energies 
x"(qu)), is mainly contributed by the asymptotics of the correlators (n (r) • rr^O)) in 
the range ^ > \x\ ^> £ s (£ t ~ v s /m t ), where /i cannot be replaced by a constant: 

(n-(r)n-(O)) ~ -L Ko (m t r), (n + (r)n + (0)> ~ -L K 2 (m t r). (17) 

Thus, even before reaching the critical point m s = 0, the 5-peak in dynamical suscep- 
tibility disappears, and the effective picture of the Haldane liquid breaks down when, 
due to softening of the singlet mode, \m s \ becomes comparable with m t . 

At m s = (|V| = 3Jj_) the singlet excitations become gapless, and the intermediate 
asymptotics ( p|) are now exact. The dynamical spin susceptibility ^L(quS) displays 
thresholds at uo = m t and u = 2m t , respectively. Near the thresholds 

- g 2 - 4m t 2 ) ^ - g 2 - m\) 

X+ ( U} *-q)~ {u2 _ q2 _± m . f/ v *m X>, ^ - q) ~ ^ _ g2 _ ^ , (18) 

The correlation function of the total dimerization follows a power law (e + (r)e + (0)) ~ 
r -1 / 4 . corresponding to the Ising order parameter at criticality. Thus the critical point 
m s = belongs to the Ising universality class (with central charge c = 1/2) and signals 
a transition to a spontaneously dimerized phase with (e + ) ^ 0. In the dimerized phase 
(m s > 0) the singlet mass gap is always the smallest indicating that the spin ladder is 
not in the regime of an effective S = 1 spin chain. 

The Haldane liquid has a nonlocal topological (string) order parameter whose 
nonzero value is associated with the breakdown of a hidden discrete (Z 2 x Z 2 ) symmetry. 
In the two-chain realization of the Haldane S=l phase the string order parameter is 
expressed in terms of the Ising order and disorder variables^]: 

lim (O string (x,y)) ~ (ax) V 2 > 2 + ^i> 2 (/i 2 > 2 . (19) 

\x-y\-*oo 

(notice that the singlet mode does not appear in this expression). Since the Ising sys- 
tems are noncritical in both phases, it follows from (|T9|) that the string order parameter 
will be nonzero in the dimerized phases as well, vanishing only at the critical points 
m t = and m s = 0. 
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Thus we have demonstrated that the Haldane spin liquid is not the only possible 
phase of a disordered magnet. The distinctive features of another - dimerized phase, 
which can be tested by inelastic neutron scattering and NM R experiments is the absence 
of coherent single-magnon modes and the 7r-periodicity of the spin excitation spectrum, 
as opposed to the undimerized phase where the ratio of the energy gaps at q ~ and 
q ~ 7i is 2. The Haldane- and dimerized phases are separated by critical points, either 
of the Ising type or belonging to the universality class of the critical S = 1 quantum 
spin chain. 

We are grateful to Yu Lu, Ilya Krive and Alexander Gogolin for valuable discussions. 
One of us (A.A.N.) acknowledges hospitality and support from the International Centre 
for Theoretical Physics. 
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